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Two-dimensional t-J model is studied by a variational Monte Carlo method, using Gutzwiller- 
Jastrow-type wave functions. Various kinds of superconducting pairing symmetries are compared 
in order to determine the phase diagram of the ground state in the full J/t-n plane. Near the 
half filling where the high temperature superconductivity is expected, the d^2_y2 wave pairing 
state is always the most stable among various symmetries. The three-site term hardly changes 
the phase diagram in this regime. In the low electron density, the extended s-type wave becomes 
a quantitatively good state for large J/t although the energy gain is small. The Gutzwiller 
wave function is shown to be the exact ground state in the low-electron-density limit for the 
supersymmetric case {J/t = 2). 
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§1. Introduction 

Symmetry of Cooper pairs is considered as a key point 
to elucidate the mechanism of high-Tc superconductivity. 
Many extensive experiments suggest that the unconven- 
tional d^2_y2 anisotropic pairing is realized in high-Tc 
oxides .HP Interesting phenomena due to this anisotropic 
pairing have been studied phenomenologically. Actually 
the two-dimensional (2D) t-J model, which is considered 
as an effective Hamiltonian for the high-Tc oxides,© will 
probably have a ground-state with the d2,2_j,2-wave pair- 
ing state for the reasonable parameters. 

The 2D t-J model is derived by regarding a singlet 
between the Cu spin and the hole's spin on the neigh- 
boring O sites (Zhang-Rice singlfit) as a mobile vacancy 
in the Heisenberg spin systcmxP Different approaches, 
such ^ jinean-field theories,oLlQLrlj'LP variational^ 
ieaij'OEB' and exact diagonalization methods Jlj'lif 
support the d2.2_j,2-wave pairing near half filling and for 
J/t - 0.4. 

In this paper we study more extensively various pairing 
states in the whole parameter space {J/t and electron 
density n) using a variational Monte Carlo method, and 
determine the phase diagram of the ground state. 

It seems natural that the strong on-site repulsion 
between up-spin and down-spin electrons favors the 
anisotropic d3.2_j^2-wave pairing. However, it is not at all 
apparent why, for example, the extended s-wave pairing 
state, which does not have the on-site pairing amplitude 
either, is not realized. 

Furthermore, many discussions on the da;2_j,2-wave 
pairing have been limited to the low doping regime (near 
the half filling) and to small values of J/t. If we consider 
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the wider range of parameter space, it is necessary to 
clarify the region where the da;2_j^2-wave is stable. For 
example, in the low-electron-density side a metallic state 
is expected, while iu ±he large- J/t region a phase sepa- 
ration is predicted.E3'EZP The relation between the phase 
separation and the d^2_j^2-wave state is not so clear. Ac- 
tually from the phase diagram of the one-dimensional 
(ID) t-J model, we expect that the^ superconducting 
state is next to the phase separation.^ But even if this 
is the case in 2D, we have to clarify the pairing symmetry 
of the neighboring superconducting state. 

Meanwhile, in the low-electron densities, there may be 
a region where the s-wave pairing becomes slahlril^ LThe 
two-electron problem is analytically solvedJlj'EflEj'LJ to 
show a stable s-wave bound state for J/t > 2. This sug- 
gests that the extended s-wave pairing state is favorable 
for J/t > 2 in the low-density region. Actually in the 
ID t-J model, such a state was expectedJlS' For finite 
electron densities in the 2D case, we have to check this 
possibility. It is far from clear how this s-wave state re- 
lates to the dx2-y2-weLve state as well as to the phase 
separation. 

For the supersymmetric case {J/t = 2), or on the edge 
of this extended s-wave regipii, we expect that the Gutz- 
willer wave function (GWF)e3' becomes the exact ground 
state for n — > 0. This fact was |predicted before in con- 
nection with the ID t-J model.H It is very interesting 
to clarify the role of the GWF in the phase diagram. 

With these aspects in mind, we focus on the following 
points in this paper. (1) How and why is the dx2-y2- 
wave pairing stabilized, but not the other pairing states 
for the high-Tc parameters? (2) The ranges of the sta- 
ble d2;2_j,2-wave state as well as other ones. Namely is 
there any possibility for the other pairing state to be re- 
alized in the phase diagram? (3) How large is the region 
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where the phase separation takes place? (4) Is there any 
region where the extended s-wave state is stable in the 
low-electron density? (5) How are the variational states 
connected to the states in the low electron density limit 
and what happens at the supersymmetric case ( J/i = 2)? 
(6) How does the three-site term affect the above results? 

In order to study these questions on the phase di- 
agram, we use the variational Monte Carlo (VMC) 
method. So far, the exact treatments like the quantum 
Monte Carlo methods have been often useless for the 
2D t'J model because of the minus sign problem. On 
the other hand, the present VMC method is free from 
such an obstacle and tr£ata_£xactly the constraint of no 
doubly occupied sites.EiEJ'E3' Accurate estimates of the 
expectation values are possible, which are indispensable 
to take advantage of merits of the variation theory. 

As variational wave functions for superconductivity, 
we use Gutzwiller-projected BCS functions which were 
originally proposed by Anderson as a resonating va- 
lence bond state. Ef It has been shown that the (1^2 _y2- 
wave pairing state is the best variational state_ai; 
Gutzwiller-projected BCS states near half filling Jlj'llil' 
In this paper, we extend previous works to include a va- 
riety of pairing symmetries in the full parameter space 
of the 2D t-J model with the three-site term. We con- 
struct a phase diagram of the ground state within the 
Gutzwiller-Jastrow-type wave functions. As will be seen 
later, the resultant 2D phase diagLam shares some fea- 
tures with the ID phase diagram.Ea' 

The mean- field theories are sometimes useful ±j(x.pr£- 
dict the behaviors of the system qualitativelyD'&ElQBlj^ 
However, they treat the constraint of no double occu- 
pancy approximately. We are going to compare the ob- 
tained phase diagram with that predicted in the mean- 
field theories with the Gutzwiller approximation. It is 
found that a qualitatively correct phase diagram is ob- 
tained, but quantitatively we need some modifications 
for the Gutzwiller approximation. 

The outline of this paper is as follows: In §2 we explain 
the model Hamiltonian and the variational wave func- 
tions. In §3 the relation between our variational states 
and the slave-boson mean-field theory is mentioned. Sec- 
tion 4 is devoted to the discussion on the metallic, mag- 
netic and superconducting states for the high-electron- 
density or low-doping regime. We study the properties 
of the stable dj:2_j^2-wave state and search for the origin 
of its stability in terms of the resonating valence bond 
theory. In §5 a phase diagram is constructed in the full 
J /t-n plane, and properties of wave functions are de- 
scribed in the low electron density. We summarize in 
§6. In Appendix A the energy of the Hubbard model is 
shown for a comparison in §4.1. The brief summary of 
the Gutzwiller approximation is given in Appendix H_ 

A part of the present results was published before.o' 

§2. Model and Trial Wave Functions 

We consider the t-J model defined in the Hilbert space 
without double occupancy of electrons. In this paper we 
include the so-called "three-site term" or "pair-hopping 
term" , Ti^ , which is usually neglected but is present 
in the effective Hamiltonian derived from the Hubbard 



model@) Thus our model is written as 
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+ _^c] ,,Cj+r',<T) , (2.4) 

where represents the sum over the nearest-neighbor 
pairs. T and r' run as vectors pointing to the nearest- 
neighbor sites. We take t as the energy unit in the fol- 
lowing. This model has four parameters: J/t, j''^) / J, 
electron density n — N /N^ and magnetization ni = 
{N^-Ni)/N^, where N (= N^+N^) is the total electron 
number, A^a the number of sites and A^^ the number of 
electrons with spin tr. We study mostly the nonmagnetic 
case {m — Q). 

The three-site term, 7^3, shows up in the effective 
Hamiltonian of the Hubbard model in the strong cou- 
pling regime {U ^ t)£B In this case J'-^-' is equal to 
At^ /U{= J). For the electron-doped high-Tc cuprate (Nd 
system), the doped electrons directly enter the Cu or- 
bitals,|-sp that the Hubbard model becomes an adequate 
modeled and thus we should assume J^^^ = J. 

On the other hand, if the effective Hamiltonian is 
derived from a hole-doped |CpP2 model including Cu 
d^2_^^nd O 2p orbitalsJa-Ej) could b e nega- 

tive.E3'E3 Therefore in this paper we use 

J(3) 

as a pa- 
rameter and discuss its effect on the phase diagram. 

For this Hamiltonian, we study Gutzwiller-Jastrow- 
type variational wave functions defined as: 



nn[i-{i-'?(^^^)}' 



(2.5) 



with Tji — \rj —r(\. The function <& is a one-body mean- 
field- type wave function; we consider (1) a simple Fermi 
sea $F as a metallic state, (2) a Hartree-Fock antiferro- 
magnetic (AF) function for the AF ordered state, and 
(3) a BCS wave function with various pairing symme- 
tries. So far these types of wave functions have been 
widely used for the small-doping region. In this paper, 
we study these states in the wholep^aiameter regime. 
The AF state state is defined asEl-ElH) 
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Here, sgn(cr) takes either +1 or —1 according as cr =1 or 
I; £k — —2t{coskx + cos ky) and the wave vector K = 
(tt, it) represents the AF reciprocal lattice vector. Aaf 
is an AF variational parameter; $af is reduced to the 
simple Fermi sea $f as Aaf — > 0. 
The BCS wave function is given aJ^^ 



with 



$sc = 



Vk_ 



n[ 



|0), 



£fc - M + Vi^k - m)^ + 



(2.7) 



(2.8) 



For most cases we take /i as the chemical potential for the 
noninteracting system, /iQ. However the optimal value of 
H ought to change due to the correlation. We will check 
the fi dependence of the variational energy in §4.3. 
The pairing symmetry of is chosen as follows: 



A, 



A 

A sin k-r 



sin fc,. 



A (cos kx 
A (cos kx 



(i) 
(ii) 
(iii) 
(iv) 



s wave 
d,.,, wave 



J-xy 

C cos ky) (iii) s+d wave 
cosfc,,) (iv) s+id wave. 

(2.9) 

In each case A (> 0) is a variational parameter and when 
A = 0, <l>sc is reduced to the simple Fermi sea, $f- Note 
here that the expectation value of the order parameter 
(cj-f c^fe|) = Asc is not necessarily equal to A. Asc is to 
be calculated as an expectation value in the variational 
state. Actually Age is largely reduced from A neaiLthe 
half filling due to the strong electron correlation. t3'EP 

By changing the value ofC (—1 < C < 1) and 9 
(0 < < tt) in the cases (iii) and (iv) , various symmetries 
are continuously realized. For example, A^ of the case 
(iii) is rewritten as: 

Afe = As(cos kx + cos ky) + Ad(cos kx — cos ky), (2.10) 

with As = A(l + C)/2 and Ad = A(l - C)/2. Thus 
the state with C = — 1 represents the pure dx2-y2 sym- 
metry, the state with C = s+d state, and the state 
with C = 1 the pure extended s state, respectively. For 
the case — 1 < C < 1, the wave function has different 
gap amplitudes in the directions of x and y axes, which 
conflicts with the symmetry of the lattice. On the other 
hand, A^ of the state (iv) represents a complex mix- 
ture of the dx^-y^ and extended s states. Especially at 
= Tr/2, Ak is the so-called s-|-«d statcEf 

Many-body effects are introduced into eq. (2.5) by a 
two-body Jastrow correlation factor 



V = 
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This factor reduces the amplitude of the wave func- 
tion according to the charge configuration. By requir- 
ing 7^(0) = 0, the local constraint of the prohibition of 
double occupancy (Gutzwiller projection) is strictly sat- 
isfied. For r 7^ 0, we consider two cases: (a) r]{r) = 1, 
which is nothing but the Gutzwiller projection Vd = 



Y[j[^ ~ "'j,T"'i-i] ^^'^ (^) intersite correlation factor 
r]{r) 
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Here L is the linear dimension of the lattice: iVa = L x 
L. We abbreviate this correlation factor as T'tll- It 
works as repulsive Jastrow correlation for > and 
attractive one for v < 0, and is reduced to Vd for — 0. 
This factor is a Tomonaga-Luttinger-liquid (TLIA type 
and has been inv£atigated in detail for the IDjOEj-* and 
the 2D systems.O Here we choose this form only for 
its simplicity as an intersite factor. Alternatively one 
can use a short-range correlation factorpj'L^ but the 
difference of the optimized variational energy between 
the two choices is negligibly small. 

§3. Relation between Variation Theory and 
Slave-Boson Mean-Field Approximation 

The phase diagram and physical properties of the 2D 
t-J model were extensively studied in the slave-boson 
mean-field approximation.B) Thus it is important to clar- 
ify the relation to the present variation theory. In the 
slave-boson mean- field approximation, the., electron op- 
erator at site i and spin a is decoupled a,a& 

Ci,a = bifida, (3.1) 

where bt is a slave-boson operator which represents the 
vacant site (holon), and fi^a- is a fermion (spinon) oper- 
ator. 

The mean-field approximation using this slave-boson 
representation gives a i-wave pairing state as the ground 
state near half filling.& ErLl^ In this state, the fermion 
(spinon) degrees of freedom have a BCS order param- 
eter 

and its wave function is the same with eq. (2.7). 

On the other hand, tke^slave bosons are assumed to 
form a bose condensate. eIlP Since the slave bosons have 
hard-core nature or two bosons cannot occupy the same 
site, it is not clear whether the bosons really form a 
bose condensate in the two (or quasi-two) dimensional 
systems. In the slave-boson mean-field approximation, 
however, the hard-core nature is neglected and thus the 
ground-state for the bosons is written as 



(&Lo)"'lO), 



(3.3) 



where iVf, is the number of holes. 

We can take account of the hard-core nature by using 
a Gutzwiller projection operator for bosons 



n(foLo)"'lo), 



(3.4) 



instead of (3.3). Even in this case, the wave-function 
amplitude for the slave boson is equal to 1 for any boson 
configurations. 

The electronic wave-function in the slave-boson 
scheme is the product of the fermion and boson wave 
functions. Therefore the ground-state obtained in the 
slave-boson mean-field theory is equal to eq. (2.7); our 
variational calculations take full account of the effect of 
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the Gutzwillcr projcatinn. 

The VMC methodcJ^EZl' enables us to estimate accu- 
rate expectation values even in the presence of the Gutz- 
willer projection. We use square lattices with L ~ 10-26 
for the metallic states and L — 10-20 for the ordered 
states; mainly the system with L = 10 is used for 
because the size dependence is negligible for L > 10 as far 
as energy is concerned. The periodic(x)-antiperiodic(?/) 
boundary conditions are used and is chosen so as to 
satisfy the closed-shell condition. To suppress statisti- 
cal errors, we take sufficient sample numbers (3 x 10^ - 
2 X 10^) and Metropolis trials (50 at maximum for each 
electron) between the samplings. 

To treat the coherent states for the superconductivity, 
two VMC schemes haveJpeen applied: one by using the 
grand canpiiical schemellj' and one by fixing the electron 
number The results of the two methods ought to 

coincide in the thermodynamic limit. The former scheme 
has a merit in calculating the order parameter directly, 
but it needs much longer CPU time. In this work we 
adopt the latter scheme. 

§4. Pairing Symmetry in Low-Doping Regime 

In this section we study the region of high electron den- 
sity or low doping, laying emphasis on the stable pairing 
symmetry of superconductivity. In §4.1 we study metal- 
lic states as a starting point, since the other states with 
AF order or superconducting order parameter reduces to 
the GWF when the order parameter vanishes. The en- 
ergy gain due to the occurrence of the order is discussed 
in comparison with the variational energy of the GWF. 
The instability of metallic states against the phase sep- 
aration and the effect of the Jastrow correlation factor 
are described. Next we investigate various pairing states 
for the half-filled case in §4.2 and for less-than-half-fiUed 
cases in §4.3. We show that the ordered states are free 
from the instability against the phase separation for a 
small value of J/t. Finally in §4.4 we look at the prop- 
erties of the stable dj.2_j^2-wave state. 

4-1 Metallic states 

First, we consider the GWF, = T^d'i'F, as a metallic 
state. In Fig. 1 three energy components: transfer energy 
Et = (Ht), exchange energy Ej = {Hj) and the three- 
site contribution E^ — (Tia) are plotted as a function of 
n; size dependence is almost negligible in this scale. 

Near the half filling, Ej is upward convex. This means 
that when J/t — s- oo the system must have a phase sep- 
aration. In order to study this kind of phase separation 
in the small- J/i region, it is necessary to study the n 
dependence of variational energies more carefully. Fig- 
ure 2 shows the close-up of this dependence near the half 
filling. On one hand, Et/t and E^/j'^'^^ are linear with 
respect to (5 = 1 — n; on the other hand, Ej / J does not 
behave linearly but 

Ej/J = Ej{S ^ 0)/J + const, x S°-\ (4.1) 

for wide range of n.il' This 5-dependence means that the 
total energy E = Et + Ej + E^ is always upward convex 
for any non-zero value of J/t. Thus the GWF is unstable 



in itself against the phase separation. 

As shown in Appendix A, this is in sharp contrast with 
the GWF in the Hubbard model, for which both Et/t 
and Eij/U are proportional to J as (5 0. The total 
energy is always downward convex, leading to no phase 
separation. This indicates that as long as the 2D GWF 
is concerned, the t-J model for J/t > is not connected 
to the Hubbard model near the half filling. 

We study this instability also for the TLL-type func- 
tioii_4'TLL = T'tll'&f, which Valcntf and Gros stud- 
ied.ES The results are summarized in Table I, indicated 
as TLL. As n — > 1 (i.e. ^ — > 0), intersite correlation 
factors generally become less effective, so that the oc- 
currence of the phase separation near the half filling 
is not affected even by this intersite correlation factor. 
Namely, as far as a form of V^f is assumed, the insta- 
bility against the phase separation is inevitable. Later 
when we study the projected BCS state, this phase sep- 
aration is removed. 

12 Half-Filled Case 

First, let us look at the results for the half filling. Fig- 
ure 3 shows the variational energy for the the case (iii) 
in eq. (2.9) for various values of C . For the d^2_j^2-wave 
case (C = —1-0), Ej/J is symmetric with respect to 
A/t = 2 and thus has two minima at A/t ^ 0.57 and 
7. (This fact, was not discussed in the previous VMC 
calculations.lU) As the parameter C increases, the min- 
imal value monotonically increases, as seen in the inset, 
where the scale is expanded. Thus the d^2„j,2 wave has 
the lowest energy in the class of variational states (iii). 
For the extended s wave (C ~ 1.0) Ihere is no depen- 
dence on A of the variational energy.EJ) This is because 
/i = at half filling and thus becomes a constant, 
which means that ^'gc is equivalent to ^'g- 

In Fig. 4 we show the variational energy for the s-|-jd- 
type wave function (iv) in eq. (2.9). As in the inset the 
energy minima for 7r/2 < < tt are degenerate within the 
error bars. This degeneracy ia due to the SU(2) symme- 
try of the Heisenberg model it was proved that under 
the Gutzwiller projection the state with Tr/2 < < tt are 
all the same state.B 

These VMC calculations for the projected BCS-type 
wave functions ^'sc are consistent with the slave-boson 
mean-field theory. However, the actual ground state at 
half filling should be an AF magnetically ordered state. 
In fact, ^I^AF — Vd'^AF is slightly stabler than ^fgc and 
the best variational state among the present states. This 
feature is i|tL.contrast with ID systems, where ordered 
states — AFE5 and s-type superconductingpl' states — 
have higher energies than the GWF. The expectation 
values of {Si ■ Sj\Sov the nearest-neighbor i and j are 
-0.3206 for *afH and -0.3199 for the A^2_y2 or thi; 
s-fid state, which are close to the best estimation,c3' 
-0.3346. 

4-.3 Symmetry of Cooper pairs in less-than-half-filling 

Now, we consider less-than-half-filled cases. We con- 
centrate on the s-f-id-type wave function, (iv) of eq. (2.9), 
which contains the dj.2_„2 and extended s-wave states. 
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All the other states give higher energy than the s+id- 
type state. For example, the s+d-type-wave states (iii) 
in eq. (2.9) with — 1 < C < 1 always give higher varia- 
tional energy as in the half-filled case. As for the d^y- 
wave state (ii), we find that each energy component is a 
monotonically increasing function of A/i for every elec- 
tron density. Thus there is no chance for the da;^ wave to 
be stabilized. The behavior of the conventional s-wave 
state (i) is very similar to the extended s wave in the 
high density regime; the energy is always higher than 
the da;2_j,2 wave for small 5. 

Let us return to the s-|-zd-type state. Figures 5(a)- 
(c) show the energy components of the s-|-id-type wave 
function. As a typical electron density, we use 5 = 0.12 
(n = 0.88). The d^2_y2 wave (9 = tt) is the most fa- 
vorable except the three-site term. Thus, when there is 
no three-site term (J*^^^ — 0), the d2.2_^2-wave state has 
always the lowest variational energy.cj' As an example, 
we show in Fig. 6 the total energy E/t for J/t = 0.25 
and several values of 9. The situation does not change 
for J/t ^ 2.0. 

Here we mention the effect of the three-site term Ha- 
The variational energy for the extended s wave is lowered 
in the presence of Tis- We have found, however, that the 
most stable state is still the (1^2 _y2 wave, even when the 
value of J(3)/J is as large as 1.0. On the contrary, the 
dx2-y2 wave is further stabilized in the ferromagnetic 
case (J^^^ < 0). The situation does not change for any 
value of J/t and for electron densities n = 0.52-0.96. 

Figure 7 shows the optimal variational energy for the 
d2,2_j^2-wave state for n ^ 1 and J/t = 0.5. In com- 
parison with the metallic GWF, there is a definite en- 
ergy gain due to the superconducting order parame- 
ter. In Fig. 7 we also plot the data of the AF state, 
^AF = 'Pd'&AF and the 7r-flux state. Just at half fill- 
ing V^AF is slightly lower than the d2.2_j,2-wave state. 
However even at (5 = 0.02 the d2.2_j^2-wave state has a 
lower energy, so that the AF state disappears immedi- 
ately when holes arcLdoped in the 2D Heisenberg system. 

The TT-fiux state,B' which is equivalent to the s-fid 
state at half filling, has exactly ^a^me energy with 



the da;2_j,2-wave state for n — In For the doped 
case, however, the vr-flux state has much higher varia- 
tional energy than the (1x2 -y2 wave. It is even higher 
than the AF state. Therefore it is not likely that the 
d2;2_j^2-wave state is mixed with the 7r-flux state. 

Finally, we refer to the dependence of the energy. 
We calculate the total energy for some J/t by varying 
the value of ^ in eq. (2.8). The smaller J/t is, the more 
weakly E/t depends on /i. As an example, the result for 
J/t = 0.5 is plotted in Fig. 8, where Aopt/i ^ 0.55. The 
energy improvement by optimal /i upon the value for /io 
is slight. Thus the use of the noninteracting value /io is 
justified for small J/t to some extent, although this is 
not the case for a value of J/t as large as 2. 

^.^ Properties of dx2_^y2 wave 

In this subsection we compare the properties of the 
d2,2_y2-wave state with other variational states in or- 
der to consider the reason of its stability. In Table I, 



we compare each energy component for some variational 
states; the optimal functions are used for the TLL, AF 
and d2.2_j,2-wave states. The differences from the GWF 
value are shown. 

First of all, we notice that the metallic state with Jas- 
trow correlation (indicated as TLL) has slightly lower ki- 
netic energy (< 1% ) than the GWF. This is because the 
Jastrow (repulsive) correlation factor tends to separate 
the electrons so as to gain the kinetic energy. However 
its gain is so small that TLL state cannot overcome the 
d:r2_j^2-wave state. 

In contrast to this, the other ordered states (AF, 
d, s-|-id state) unanimously gain the exchange energy, 
namely {Si ■ Sj) and {—niUj/A), but all of them lose 
the kinetic energy. As the doping increases the energy 
gain for the exchange term becomes smaller and smaller. 
Finally at some doping, the energy gain due to the AF 
order or superconducting order vanishes. This tendency 
gives the n-dependences of the variational energy shown 
in Fig. 7. 

Comparing the d2:2_j^2-wave with the AF states, we 
can see that the loss of the kinetic energy is always 
smaller than that for the AF state, while the gain of 
the exchange energy is larger for the d^2_j,2 wave for 
every electron density. This makes the dx2^y2 the best 
variational state. 

In Table I, we also include the energy components 
for the s+id and extended s states. For comparison 
we borrow the same value of A/t with the optimal 
d,j.2-y2 case, namely A/i = 0.7,0.55,0.40, and 0.25 for 
n = 096,0.88,0.80, and 0.72, respectively. As for the ex- 
change term {Si ■ Sj), the da;2_j,2-wave state has a large 
energy gain. For example at n = 0.88 {5 — 0.12) it gains 
29% of the exchange term of the GWF. The s+id- wave 
also gains the exchange term similar to the d^2_j^2 wave, 
but its magnitude is slightly smaller. For the extended 
s wave the energy gain is very small. This is the reason 
why the d2,2_j^2-wave state is favored by the exchange 
term J. We discuss this phenomena in terms of the res- 
onating valence bond (RVB) representation of the wave 
functions. 

Equation (2.7) can be rewritten as 



^'sc = j]^ ■ Gxp 



|0), 



(4.2) 



where a^j- is Fourier transform of Ofc. If we fix the elec- 
tron number by projecting out the N electron states, we 
obtain 



k 



1 



Et t 



|0). (4.3) 



This form of the wave function is the superposition of 
various RVB configurations. Each valence bond (VB) is 
represented by cJ|C^|, whose amplitude is a^.j. 

By investigating the actual form of a^j from the def- 
inition in eq. (2.8) for various symmetries, we can see 
the difference between the dx2-y2 wave and s-|-id, or ex- 
tended s wave. The main difference is that, among the 
dominant a^j 's, the next-nearest neighbor bond exactly 
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vanishes for the dx2_y2 state, i.e. 



A, 



= 0, 



£fc - M + Vi^k - f^V + 

(4.4) 



because of the antisymmetric nature for exchanging 
and ky. On the other hand, for the s+id and extended 
s cases, ai^i^x+y has a comparable amplitude with the 
nearest-neighbor ai^i+x or ai^i+y. 

Due to this main difference, the RVB amplitude for 
the nearest- neighbor bond for the dx^-y^ state is more 
enhanced than the s-Kid or extended s state. As a re- 
sult, the d2;2_j,2-wave state is the most favorable for the 
exchange energy. (Reversely we can say that the next- 
nearest-neighbor exchange energy favors the extended s- 
wave state.) 

To see this more closely, the spin structure factor. 



iq{rj -re) 



(4.5) 



is compared in Fig. 9 for the optimized dx2_y2 wave and 
extended s wave. The antiferromagnetic correlation is 
enhanced for the dx2-y2 wave, which leads to the gain in 
the exchange energy. 

Next we consider the difference in the kinetic energy in 
Table I. The ordered states always lose the kinetic energy 
compared with the GWF. Without taking account of the 
Gutzwiller projection, we can estimate 



1 £fc ~ M 

k 



{cos kx + COS ky) . (4.6) 



In the case of A = 0, (ej, — fi)/Ek = sign(efc — fj,), but in 
the presence of A, the amplitude of (c||Cj,|) generally 
reduces. This causes the loss of the kinetic energy for the 
ordered state. Since there is a zero point at kx = iky 
on the Fermi surface for the d2:2_j,2-wave state, the loss 
is the least among the superconducting states. This fact 
also favors the dx2-y2 symmetry. 

Finally we discuss the {—riinj /A) term. For example, 
the GWF gives {-n^^lj/^) = -0.1925 for n = 0.88 case. 
This is understood from the fact that {rii) — 0.88, then 
(nj)2/4 = 0.1936. As shown in Table I, the value of 
this term for the other states does not change so much 
from that of the GWF. For example, {—fiirij/A) changes 
only 0.5% for the dx2_y2 wave. Actually the charge- 
density correlation function in the real space shows that, 
in the dx2-y2 case, the amplitude of the nearest and 
next-nearest neighbors is slightly enhanced. This is in 
contrast to the GWF which shows dips at those sites 
due to the exchange hole. Meanwhile, the extended s 
wave scarcely modifies the GWF value. Thus, we expect 
that an attractive correlation factor gives more energy 
gain. We check the effect of the intersite correlation fac- 
tor by introducing eq. (2.11) to the da;2_j,2-wave state: 
'?'tll^SC- However, we have found that the energy is 
hardly affected by the intersite correlation. Therefore 
the stability of the dx2_y2 state is mainly caused by the 



energy gain of the exchange term which is understood 
from the RVB picture as discussed above. 

§5. Phase Diagram 

In this section, we construct a phase diagram in the 
whole J /t-n plane. First of all we show the summary 
of the phase diagram in Fig. 10. In the following sub- 
sections we describe the details of determining the phase 
boundaries. 

In §5.1 we discuss the phase boundary between the 
d^2_y2-wave state and the metallic states. The stabiliza- 
tion of the extended s wave in the low electron density 
is studied in §5.2. The properties of the GWF around 
the supersymmetric case ( J/i = 2) is discussed in §5.3. 
In §5.4 the phase separation in the large J/t region is 
studied. We consider in §5.5 the effect of the three-site 
term on each phase. 

5.1 Range of dx2_y2 wave 

To begin with, we determine the phase boundary be- 
tween the d2,2_j^2-wave state and the metallic state. Here 
we use the GWF as the metallic state, because the en- 
ergy gain due to the intersite Jastrow correlation factors 
is very small. Figures 11 and 12 show the energy com- 
ponents of the da;2_j,2-wave function for various electron 
densities ranging from n = 1 down to 0.04. 

The kinetic energy Et/t shown in Fig. 11 goes to 
as the electron density approaches either or 1. The 
latter is due to the Gutzwiller projection, Vd- Further- 
more Ef/t is a monotonically increasing function of A/t 
in every density. This means that the kinetic energy 
does not favor the superconductivity for every electron 
density. On the other hand, Ej/J shown in Fig. 12 de- 
creases at finite values of A/t. The exchange term gives 
the attractive force in the dx2^y2 channel as expected. 

Thus the dx2-y2 state is not realized without the ex- 
change interaction {J/t = 0). For finite J/t, the varia- 
tional energy 



E = tx {Et/t) + J X {Ej/J) 



(5.1) 



has a minimum at a finite value of A/t. The obtained 
optimal values of A/t are plotted in Fig. 13. From this 
figure the phase boundary is estimated for every elec- 
tron density. For example, 5 = 0.2 {n — 0.8) case the 
phase boundary is at J/t = 0.1. This boundary is shown 
by a solid line with solid triangles in Fig. 10. For the 
higher doping case n ^ 0.64 [Fig. 13(b)], it is difficult 
to determine the boundary definitely, because small A 
seems to be realized for small value of J/t. However, this 
residual small A is not so meaningful practically, because 
the energy gain is as small as that of the TLL (metal- 
lic) state. Therefore, we estimate the phase boundary 
by extrapolating the straight fitting lines as indicated in 
Fig. 13(b). The boundary thus determined is represented 
by a dashed line with reversed triangles in Fig. 10. 



5.2 Range of extended s wave 

Next we discuss the stability of the extended s wave 
in the low density region. As-shDWii-hy the analytic 
treatments in the n ^ Q hmit,E£l'E3'E3'EI the free elec- 
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tron state is unstable against the s-wave bound state for 
2 < J/t < 3.4367. We will discuss this possibility of the 
extended s wave in the finite electron densities. 

In Figs. 14 and 15 we show the expectation values of 
energy components for the extended s wave. [The result 
of s-wave symmetry (i) of eq. (2.9) is quantitatively sim- 
ilar to that of the extended s wave.] The kinetic energy, 
Et/t, rapidly approaches zero as A/i increases. As re- 
gards Ej/ J ^ energy gain is hardly seen for small doping 
as shown in Fig. 15(b). This is consistent with the result 
in the half filling where Et/t does not depend on A as 
discussed before. On the other hand for small electron 
densities [Fig. 15(a)], Ej/ J appreciably decreases. This 
is the cause of the stability of the extended s wave in the 
low density. 

In Fig. 16 we plot the total energy for n = 0.04 as an 
example. No energy gain from the GWF can be seen for 
J/t < 2, while the energy gradually decreases for J/t > 
2. This is consistent with the analysis of two electron 
systems, where the electrons form a s-wave bound state 
for J/t>2. 

For finite electron densities, the stable region of the 
extended s wave is shown in Fig. 10. Since the energy 
gain is very small and the same order with that obtained 
by a metallic state with an intersite correlation factor, it 
is not easy to determine accurately the boundary of the 
extended s wave. 

In Fig. 16 we also plot the energy gain due to the 
intersite correlation factor, ^'tll- For J/t < 2.0, the 
optimal value of i' in eq. (2.12) is positive, which cor- 
responds to the repulsive Jastrow factor. For example, 
we get i/opt ~ 0.25 (0.20, 0.14, 0.06) for J/t = (0.5, 
1.0, 1.5). Just at J/t = 2.0 we find that the Jastrow 
factor does not reduce the variational energy, which in- 
dicates that the GWF is a good variational function for 
this "supersymmetric" case as in ID.cJ We will come 
back to this point shortly. For J/t > 2.0, I'opt is nega- 
tive, namely, t^opt 0.09 (-0.19, -0.27) for J/t = 2.5 

(3.0,3.5). This means that jfiffective attractive interac- 
tion exists between electron.O 

5.3 Supersymmetric Case 

Now we investigate the supersymmetric case {J/t = 2) 
more in detail. For this special parameter, ±he GWF is 
shown to be exact in the low density limit .£3' Actually, 
the two-electron problem can be solved analytically. We 
consider the GWF with two electrons. 



(5.2) 



If H (including H3) is applied to this state, the resultant 
equation is written as. 



-2zf^G + {2t-J-{z-l)J^^^E, (5.3) 



where z is the number of nearest neighbor sites (here 
z = 4) and 



This means that if the condition 

2t- J-{z- l)j(^) = 



I.e. 



J 
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(5.4) 



(5.5) 



(5.6) 



1 + (Z-1)J(3)/J 

is satisfied, the GWF is an exact eigenstate with energy 
-2zt, irrespective of A^a [ in fact (*g|S)/(*g|*g) = 
0{1/Nii) ]. The energy is identical with the noninteract- 
ing case. As will be mentioned in §5.5, the free electron 
state has no instability against the bound-pair states of 
both s and d waves when eq. (5.6) holds. Thus the GWF 
is shown to be the exact ground state. This property 
does not depend on the lattice form and dimension. 

One carLinterpret this as the same aspect with the ID 
t-J modeLo-* Namely the electron hopping term, which 
favors the electrons mutually apart, and the exchange 
term, which works attractively, are well balanced and a 
kind of "noninteracting" state is realized. 

5.4 Region of phas^ separation 

As was shownjlj' ili a phase separation appears in the 
2D t-J model for the large J/t region. Here we discuss 
it from our variational calculations. The phase sepa- 
ration occurs between the Heisenberg system (without 
hole) and a hole-rich phase. For the small-n region, the 
state with phase separation is approximated as a state 
with the Heisenberg island and the remaining vacuum. 
The energy of this state is 



E 



PS 



uEta 



-1.1693nJ, 



(5.7) 



where £^Hcis.ri^ the ground state energy of the 2D Heisen- 
berg modeI.ES Comparing this energy with the opti- 
mized variational energies of the GWF, the d^2_y2-wave 
and extended s-wave states, we determine the phase 
boundary, which is also shown in Fig. 10. 

In Fig. 16, the energy Epg is also plotted; since the 
J/t-dependence is strong, we can determine the phase 
boundary accurately in the small electron density region. 
Table II summarizes the phase boundary (Jc/t) obtained 
in this way. The extrapolation of these values to n = 
leads us to Jn/± = 3.437, which is quite close to the exact 
value 3.4367M 

As shown in Fig. 16 the GWF is a good trial function 
in the low-electron density and the variational energy 
has weak J/t dependence. Therefore we can broadly esti- 
mate Jc/t a,s n ^ from the GWF. In the limit of n ^ 
the energy components of the GWF are analytically ex- 
pressed as Et/t = -An, Ej/J = -n^, Ej./j'^^'^ = -in^. 
These values are the same with the free electron sys- 
tem, because double occupation can be substantially ne- 
glected. Using -Egwf — and _Eps in eq. (5.7) we 
obtain the phase boundaryta* 



Jc/t-4/|£;Heis./J| =3.42. 



(5.8) 
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As n increases, our estimation of the phase bound- 
ary becomes inaccurate. For one thing, we have as- 
sumed that the phase separation takes place between the 
Heisenberg island and the vacuum. Another thing is that 
we have compared variational energies with the accurate 
energy of the phase separation with i?Heis- Nevertheless 
our estimatioD of the phase boundary agrees well with 
other resultstlP up to fairly high density {n ^ 0.9). 

5.5 Effect of the three-site term on the phase diagram 
Next, we consider the effect of the three-site term 
on the variational energies and the phase diagram. In 
Figs. 17 and 18, i^a/j'^' is shown for the d2.2_j^2-wave 
and extended s-wave states, respectively. For the Ax2-y2 
case, i^a/J^^-' behaves similarly to Et/t, while for the ex- 
tended s wave it decreases with increasing A/t and its 
shape resembles that of Ej / J rather than Ei/t. 

First we consider the the boundary between the (1^2 _y2 
wave and the metallic phase. In Fig. 19 we plot the 
Aopt/i of the Ax2-y2 wave for some values of J^'^'^ / J 
and n — 0.8. The phase boundary of the (lx'^_y2 wave, 
Jc/t, where Aopt vanishes, scarcely changes from that 
for J(3) = 0. This aspect is common for other densities 
{n = 0.52-0.96). Thus, the region of the (Ix2-y2 wave 
hardly changes at least for n > 0.64 by the three-site 
term. 

The phase boundary of the extended s wave, on the 
contrary, depends on the three-site term considerably. 
For the extended s wave, E^/J^^^^ decreases with increas- 
ing A/t as shown in Fig. 18. Thus the extended s wave 
becomes stabler for J^^^/J > and its region expands. 
The shift of the boundary is also shown in Fig. 10 by 
an arrow for Es/J'^^^ = 0.5. Inversely, for j(3)/J < 
the energy of the extended s wave increases so that its 
region disappears for the case of j'^'^^ / J — —0.5, namely 
the GWF is stable for every value of J/t. 

This behavior is consistent with the analysis of two- 
electron problem. If we solve the two electron problem 
including Ti.^, we obtain the eigenvalue equation 



-^(-)=^(1 + 3J(^)/J) 
na \a J 4r 













- 1 


TT \ 







(5.9) 



for the twOj-electron bound state with the extended s 
symmetry,c3' where a — —E/At > 1 and K{k) is the 
complete elliptic integral of the first kind. Equation (5.9) 
has a solution when 



i(l + 3j'»'/^)>i 



(5,10) 



As a result, for fixed J^'^^'/J, the critical value above 
which the two electrons form a bound state is 

Jc 2 , ^ 

(5.11) 



t 



1 



This value is exactly the same with the condition eq. (5.6) 
where the GWF becomes an exact eigenstate. For 
J(3)/J = 0.5, this critical value is J/t = 4/5. This 
value coincides with the boundary shown in Fig. 10 
which we determine from the VMC calculations. When 
j(3)/J < —1/3, there is no solution for the two-electron 
bound state. 



The phase boundary to the phase separation is also 
estimated. The shifted boundaries for J^^-'/J — ±0.5 are 
shown with open circle and dotted line in Fig. 10. For 
j(3)/j = 0.5 the extended s-wave state becomes so stable 
that the phase separation does not take place for small 
n. On the other hand, for J*^^) J j — —0.5 the extended s- 
wavc state is not stabilized and thus the phase separation 
expands to smaller J/t region. 

§6. Summary and Discussions 

By using the variational Monte Carlo method, we have 
studied the t-J model with the three-site term for the 
square lattice, and obtained a phase diagram in the 
full J/t-n plane. We have used Gutzwiller-Jastrow-type 
wave functions describing metallic, antiferromagnetic, 
and various singlet-pairing superconducting states: the 
BCS-like s, (lxy^ s-|-d-type, and s-|-id-type waves, which 
include the (1^2 _y2 and extended s waves. 

Main results are summarized in the phase diagram 
shown in Fig. 10. Remarkable points which we have con- 
firmed and found are: 

(1) At half filling, many pairing states including the 
dx2-y2 and the s+id-wave states are degenerate and the 
most stable, in accordance with the analysis of the lo- 
cal SU(2) symmetry.cll' B Nevertheless, the antiferromag- 
netic state is more stable here. 

(2) For plausible parameters of high-Tc superconduc- 
tors, the d2.2_j,2 wave is the most stable state among var- 
ious pairing symmetries. The range of the da.2_y2 wave 
hardly changes even if the three-site term is added. Sub- 
stantially, there is no chance to realize s+d, s+id, and 
extended s waves in this region. 

(3) The dxy wave always has higher energy than the 
GWF except for extraordinarily large negative value of 

(4) In the region of low electron density and J/t > 2.0, 
the extended s wave becomes a good state, which almost 
reproduces the exact value of the phase boundary for 
n ~* 0. However the energy gain from the GWF is small. 
This phase corresponds to the low-electron-density phase 
predicted in the ID t-J model, which consists of two- 
electron singlet bound-pairs. 

(5) Metallic states V^f are unstable in itself against 
the phase separation for n ^ 1 and J/t > 0. This is in 
contrast with the cases of the ID t-J, and ID and 2D 
Hubbard models. 

(6) For the supersymmetric case {J/t = 2), the GWF 
is exact in the n — > limit and behaves like the non- 
interacting state in th£_.low-electron-density region just 
hke the ID t-J modelB) The GWF is also a very good 
trial wave function in the region up to n ^ 0.25. 

Our result that the dj.2_y2 wave is stable for the 
parametcjp^ of. high-Tc oxides confirms the early VMC 
results.Eij'Ell'H The obtained phase diagram is consis- 
tent with the one by the high-temperature expansion,0^ 
and is qualitatively similar to that of the VMC calcu- 
lations by Dagotto et altB We suspect, however, that 
the qualitative and/or quantitative differences between 
the previous VMC worksOcJ'EJ' and the present one 
are attributed to some insufficient or erroneous numeri- 
cal treatments of the the former works. 
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We should not consider the results (l)-(5) are com- 
pletely decisive because of the approximate trial func- 
tions used. However, the results of the relative stability 
among various pairing symmetries are persuasive, be- 
cause we do not give partial treatment to a specific sym- 
metry. 

We also calculate the variational energies using the 
Gutzwiller approximation. EP In contrast to the VMC 
method, the Gutzwiller approximation assumes the ef- 
fect of the projection as some statistical weights. The 
details are summarized in Appendix B. Figure 20 shows 
the obtained variational energy for various superconduct- 
ing states including the d^2_y2 and extended s waves for 
n ~ 1 and J/t = 0.5. The variational energies for the 
metallic GWF is also plotted. This figure is to be com- 
pared with Fig. 7 where the same calculation is carried 
out using the VMC method. The Gutzwiller approxima- 
tion gives a similar result. In Fig. 20 there is a small re- 
gion in the vicinity of rt = 1 where s+id state is realized, 
which was found before. LP This is due to the violation of 
the SU(2) symmetry in the approximation. 

Figure 21 shows the phase diagram determined from 
the Gutzwiller approximation. As expected from the 
good coincidence of the variational energy with the VMC 
results, the obtained phase diagram resembles Fig. 10 de- 
termined from the VMC method. Thus, the Gutzwiller 
approximation gives broadly reasonable results. 

In the following we compare our results with other 
works. 

As for the stability of the dj,2_j,2-wave pairing, our 
result is common with one for the 2D d-p model. Us- 
ing similar variatitmal wave functions for this model, 
Asahata and Oguricj^ recently showed that the dx2_y2- 
wave pairing state appears in the borders of metallic and 
AF regions in the phase diagram. On the other hand, 
concerning the possibility of the extended s wave, it is 
pointed out that modification of the Fermi surface struc- 
ture in $ of eq. (2.5) is important. 

As shown by previous VMC calculations ,£3 the ferro- 
magnetic region with full momentum corresponding to 
Nagaoka's theoremO exists in the region 0.6 ^ n < 1 
and J/t ^ 0.1. On the other hand, according to the 
high-temperature-expansionj^ partial ferromagnetism 
appears in the similar region. In this paper, we have not 
discussed this ferromagnetism in detail, since it is con- 
fined in a small parameter range. The relation among 
the above two calculations and our present results is an 
interesting issue lefLfor future studies. 

Kagan and RiceH showed that the p-wave pairing is 
the leading instability for small values of J/t and low- 
electron densities by means of the i-matrix technique. In 
this paper we have not discussed this jiassibility, because 
its energy gain seems extremely smallO^ It is rather dif- 
ficult to detect a tiny energy gain by the ordinary VMC 
calculations; moreover such a small energy gain is easily 
overcome by introducing the Jastrow correlation factors 
in the metallic states. 

Another interesting issue is the possibility of spin- 
charge separation in the 2D systenis. From the density 
correlation function, Putikka et alE3 discussed the spin- 
charge separation. However as claimed later, it is not 
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yet conclusive.!^) To study TLL states for the metallic 
phase, a VMCLpiethocl21r and a Green function Monte 
Carlo methodE^P have been also used. In this paper we 
have taken account of this possibility by adopting a TLL- 
type correlation factor. Consequently, the energy gain 
thereof is very small and it is not likely that the ground- 
state becomes a Tomonaga-Luttinger liquid state in the 
present variational scheme. 
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Appendix A: The GWF for the Hubbard Model 

In this appendix, we show the behavior ol E/t as a 
function of n (n ~ 1) for the 2D Hubbard model: 

Hhub = -t {4.aCj,<y + S>Ci,<T) +^Y1 ^n^jl^ 

(A-1) 

with the GWF: 

*g(5) = n [1 - (1 - 9)nnn,^] c&p. (A-2) 
j 

The optimized results near the half filling is shown in 
Fig. 22; E/t behaves linearly as 5 ^ 0, although size 
dependence becomes a little larger as n approaches 1. 
This linear behavior is in sharp contrast with the result 
of the t-J model as discussed in §4.1. 

This discordance may be intuitively understood as fol- 
lows. By neglecting higher-order terms, H of eq. (2.1) is 
a canonical transformation of the Hubbard Hamiltonian: 

H = e-'^^Hhub e'""- (A-3) 

Therefore, the expectation value with respect to the 
GWF is written as 

= (^iWhubl^-) (A-4) 

where ^ = c^'^^'q. Consequently, we should adopt mod- 
ified functions ^I* for the Hubbard model to obtain a sim- 
ilar behavior. 

Nevertheless, iiJ-ip, where the modified wave functions 
are much better,L3'Ef|) there exists no conflict between 
the two models even within the GWF. And this GWF- 
result is consistent with the exact analysis,E3' which 
shows that the Hubbard model is connected to the t-J 
model for every n in the strong coupling limit. Anyway, 
it is necessary to confirm by some reliable method that 
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the t-J model is connected to the Hubbard model also 
for 2D metallic states. 

Appendix B: Gutzwiller Approximation 

In this Appendix wc summarize the Gutzwiller approx- 
imation for the 2D t-J model. This approximation is 
extensi|ZEly studied in connection with the VMC calcu- 
lations .eP In order to compare with the present results, 
we only have to include the estimation of the {niUj/A) 
term. 

In the Gutzwiller approximation, the expectation 
value 



is approximated b 



25 



(^o|4<tCj,<t|V'o), 



(B-1) 



(B-2) 



1 + S 

where S is the doping rate {6 = 1 — n). Similarly the 
exchange term is given by 

4 



:{^o\S^■SJ\^o)■ (B-3) 



(1-f <5)2 

For the {niTij/A) term, we introduce the following ap- 
proximation, 



i'ipol'Pdn^nj'PdliJ-'Q) 



{1 + sy 



(1+^)2 

(V'o|«i"j|V'o) 



+ (1 - nii)n^^nji{l - n^i) 
+ (1 - rij|)nj|nj|(l - Uji) 
+ (1 - nii)n^inji{l - njt)IV'o) 

(1 - - r^x)(V'o|?^^f^J#o) 

(B-4) 



Using these approximations, self-consistent equations 
are given by 

A. =^5:7.-.':^, (B-5) 



2Ek' 



6 =— — 



Ik-k' 



(B-6) 
(B-7) 



where 



7fc= 2(cos fcj, -I- cos ky), 
26 4 
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1 + S' -"^ (l + <5)2' - ■ Na d6 

These self-consistent equations are numerically solved for 
the da;2_|_y2, s-|-id and extended s-wave states. 
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Fig. 1. Energy components of the Gutzwiller wave function as a 
function of electron density. Each Symbol represents a system 
size: solid circle {L = 26), upward triangle (20), downward tri- 
angle (16), circle (14), square (12) and diamond (10). Curves 
are a guide for the eyes. 



Fig. 2. Data of Fig. 1 are re- plotted. Et/t and E^/J^^^ are plot- 
ted versus 1—n in (a), and Ej/J versus (l-n)0-7 in (b). Symbols 
have the same meanings with Fig. 1. 



Fig. 3. Energy expectation values at half filling for the case (iii) 
(s4-d-type wave) of eq. (2.9) for several values of C (—1.0, —0.6, 
—0.4, 0, 0.4, 0.6, 1.0). The inset is the magnification of the 
minimum area around A/t = 0.57. Symbols for the value of C 
are common in the inset. The system used is L = 10, and the 
sample number is 3 - 5 X 10^. Thepwilue of Green function Monte 
Carlo method is E/J = -1.1693.E|) 



Fig. 4. Energy expectation values at half filling for the case (iv) 
(s-fid-type wave) for several values of 9 (0, vr/12, 7r/4, 71/2, 37r/4, 
7r). The inset is the magnification of the minimum area around 
0.5 - 2. Symbols of 9 are common with the inset. The system 
used is L = 10, and the sample number is 3 - 5 X 10*. 



Fig. 5. Expectation values of three energy components for s-l-id- 
type wave and for some value of d (0, 7r/4, 7r/2, 37r/4, tt) in the 
same energy scale. Electron density is chosen at n = 0.88. The 
system used is L = 10, and the sample number is 3 - 5 X 10*. 
Since there are large statistical fiuctuations in Ej / J and E^/J^^^ 
for small 9 and large A/t, we show the data of two Monte Carlo 
runs for 9 = 0. 



Fig. 6. Total energy of s-|-jd-type wave for some value of 9 (0, 
7r/4, it/2, 3n/4, tt). J^^)/ J = 0. The inset is the magnified 
figure of the minimum area. Symbols of 9 are common with the 
inset. The system used is L = 10, and the sample number is 3 - 
5 X 10* 



Fig. 7. Comparison of total energies among four variational 
states as a function of n near the half filling (j(3) = 0). Symbols 
have the same meaning with Figs. 1 and 2. The arrow on the 
right axis represents the accuifate value obtained by the Green 
function Monte Carlo method.Ej' 
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Fig. 8. Total energy of the d^2_y2-wave state versus /i/t for some 
values of A/t. Note that the scale is small. The arrow on the 
horizontal axis indicates the noninteracting value of ^. 3 X 10'' 
samples are used. L = 10. 



Fig. 16. Comparison among some variational energies as a func- 
tion of J/t for n = 0.04 for the t-.J model ( jt^) = 0). For J/t < 2 
the optimal extended s wave is reduced to the GWF. Note that 
the scale is small. L = 10 and the sample number is 2 X 10^. 



Fig. 9. Comparison of spin structure factors between the GWF, 
the optimized d^2_y2 wave (A/t = 0.55), and the extended s 
wave (A/t = 0.55). The path of q is r(0,0) X{tt,0) 
M(tt, tt) F. The system with L = 10 and 5 X 10^ samples are 
used. 



Fig. 17. Expectation values of three-site term of the d^2_y2 wave 
for (a) low density (n = 0.04 - 0.52) and (b) high density {n = 
0.60 - 0.96). L = 10 and the sample number is 3 - 5 X 10*. 



Fig. 10. Phase diagram for the ground state of the 2D t-J-type 
model. Basically it represents the case of Jt^*) /J = 0. The shifts 
due to finite /J are shown b-Mhadowed arrows with their 
values. The exact phase boundarytlP between the s-wave bound 
state and the phase separation for n — » is shown by a solid 
arrow on the abscissa. For detailed explanations, see the text. 



Fig. 18. Expectation values of the three-site term of the extended 
s wave for (a) low density (n = 0.04 - 0.52) and (b) high density 
(n = 0.60 - 0.96). L = 10 and the sample number is 3 - 5 X 
10*. For the large value of A/f large statistical fluctuations are 
observed, and the results of two trials are depicted for some 
values of n. 



Fig. 11. Expectation values of transfer energy of the d^2_y2 wave 
for (a) low density (n = 0.04, 0.12, 0.16, 0.24, 0.32, 0.36, and 
0.44) and (b) high density (n = 0.52, 0.60, 0.64, 0.72, 0.80, 0.88, 
and 0.96). L = 10 and the sample number is 3 - 5 X 10*. 



Fig. 19. Optimized value of A/t for some values of J^^^ / J as a 
function of J/t. Solid lines are guides for eyes. The wave function 
has the d^2_j^2-wave pairing symmetry. 



Fig. 12. Expectation values of exchange energy of the d^2_y2 
wave for (a) low density (n = 0.04 - 0.52) and (b) high density 
(n = 0.60 - 1.0). L = 10. The sample number is 3 - 5 X 10*. 



Fig. 20. Comparison of total energies among three variational 
states obtained in the Gutzwiller approximation as a function 
of n near the half fiUing = 0). 



Fig. 21. Phase diagram obtained in the Gutzwiller approxima- 

, , r A ; r , I r tlou for thc grouud state of the 2D t- J-type model. 

Fig. 13. Optimized value of A/t for the a^2_y2 wave as a func- 
tion of J/t (a) in the high density (n = 0.64 - 0.96) and (b) in 
the medium density (n = 0.32 - 0.60). In (a) we fit the data with 
curves, and in (b) with straight lines. The arrow on the vertical 
axis in (a) is the value for n = 1.0. 

Fig. 22. Total variational energy for 2D Hubbard model as a 
function of electron density near the half filling, with respect 
to the Gutzwiller wave function. Used lattice sizes are indicated 
by the same symbols with Fig. 1. Lines are a guide for the eyes. 
3 X 10* samples are used. 

Fig. 14. Expectation values of transfer energy of the extended s 
wave for (a) low density (n = 0.04 - 0.44) and (b) high density 
(n = 0.52 - 0.96). L = 10 and the sample number is 3 - 5 X 10*. 



Fig. 15. Expectation values of exchange energy of the extended 
s wave for (a) low density (n = 0.04 - 0.52) and (b) high density 
(n = 0.60 - 1.0). L = 10 and the sample number is 3 - 5 X 10*. 
Since for large A/t statistical fluctuations are large, the results 
of two trials are depicted for a few values of n. 
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Table I. Differences of energy components from the GWF for 
J/t = 0.5 in the high electron density. The values of the GWF 
are given as the references. The system with L = 10 is used. 



The last digit for each value includes 


some error. 




n 


State 


Et/t 




riiUj 
4~ 


Total 




GWF 


0.0 


-0.2706 


-0.25 


-0.5206 




TLL 


0.0 


0.0 


0.0 


0.0 


1.00 


AF 


0.0 


-0.0500 


0.0 


-0.0500 




d^2_y2 


0.0 


-0.0493 


0.0 


-0.0493 




Exact 


0.0 


—0.0641 


0.0 


—0.0641 




GWF 


-0.1080 


-0.2216 




n ^^Q7 
— u.ooy 1 




TLL 


~ 


~ 


~ 


-0.0001 




AF 


+0.0098 


-0.0543 


-0.0000 


—0.0447 


0.96 


d^2_j,2 


+0.0053 


-0.0603 


-0.0004 


-0.0555 




s+id 


+0.0075 


-0.0568 


-0.0004 


-0.0497 




Ext.s 


+0.00vi7 


+0.0011 


-0.0000 


+0.0048 




GWF 


-0.3089 


-0.1746 


n 1 Q9^ 


— U.U 1 oo 




TLL 


-0.0004 


+0.0001 


+0.0002 


-0.0005 




AF 


+0.0173 


-0.0268 


-0.0000 


-0.0100 


0.88 


d^2_y2 


+0.0149 


-0.0499 


-0.0012 


-0.0367 




s+id 


+0.0203 


-0.0435 


-0.0011 


-0.0242 




Ext.s 


+0.0212 


+0.0007 


-0.0003 


+0.0216 




GWF 


-0.4834 


-0.1426 


-0.1575 


-0.7834 




TLL 


-0.0022 


+0.0005 


+0.0005 


-0.0012 




AF 


0.0 


0.0 


0.0 


0.0 


0.80 


d^2_j,2 


+0.0195 


-0.0357 


-0.0017 


-0.0179 




s+id 


+0.0232 


-0.0273 


-0.0014 


-0.0055 




Ext.s 


+0.0297 


-0.0005 


-0.0005 


+0.0288 




GWF 


-0.6290 


-0.1193 


-0.1257 


-0.8743 




TLL 


-0.0042 


+0.0009 


+0.0010 


-0.0021 




AF 


0.0 


0.0 


0.0 


0.0 


0.72 


d^2_j,2 


+0.0102 


-0.0143 


-0.0010 


-0.0050 




s+id 


+0.0168 


-0.0118 


-0.0009 


+0.0040 




Ext.s 


+0.0226 


-0.0002 


-0.0004 


+0.0219 



Table IL Critical values of J/t of the phase boundary between 
homogeneous and separated phases for the low electron density. 
The last digit for each value of Jc/t includes some errors. We 
represent data for L > 14. j'^' /J i=nO. The exact value for 
n = is by Hellberg and Manousakis.E3) 



n 


N/N^ 


Jc/t 


State 


0.0 


2/oo 


3.4367 


Exact 


0.0059 


4/676 


3.432 


Ext.s 


0.01 


4/400 


3.429 


Ext.s 


0.0156 


4/256 


3.423 


Ext.s 


0.0204 


4/196 


3.418 


Ext.s 


0.03 


12/400 


3.389 


Ext.s 


0.0612 


12/196 


3.345 


TLL 


0.0816 


16/196 


3.313 


TLL 


0.0938 


24/256 


3.284 


TLL 



